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Using analytical formulas as well as a finite-difference scheme, we investigate the magnetic field 
dependence of the energy spectra and magnetic edge states of HgTe/CdTe-based quantum wells 
in the presence of perpendicular magnetic fields and hard walls, for the band-structure parameters 
corresponding to the normal and inverted regimes. Whereas one cannot find counterpropagat- 
ing, spin-polarized states in the normal regime, below the crossover point between the uppermost 
(electron-like) valence and lowest (hole-like) conduction Landau levels, one can still observe such 
states at finite magnetic fields in the inverted regime, although these states are no longer protected 
by time-reversal symmetry. Furthermore, the bulk magnetization and susceptibility in HgTe quan- 
tum wells are studied, in particular their dependence on the magnetic field, chemical potential, and 
carrier densities. We find that for fixed chemical potentials as well as for fixed carrier densities, 
the magnetization and magnetic susceptibility in both the normal and the inverted regimes exhibit 
de Haas-van Alphen oscillations, whose amplitude decreases with increasing temperature. More- 
over, if the band structure is inverted, the ground-state magnetization (and consequently also the 
ground-state susceptibility) is discontinuous at the crossover point between the uppermost valence 
and lowest conduction Landau levels. At finite temperatures and/or doping, this discontinuity is 
canceled by the contribution from the electrons and holes and the total magnetization and suscep- 
tibility are continuous. In the normal regime, this discontinuity of the ground-state magnetization 
does not arise and the magnetization is continuous for zero as well as finite temperatures. 

PACS numbers: 73.63.Hs,73.43.-f,85.75.-d 
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I. INTRODUCTION 

In recent years, much attention has been devoted to 
the field of topological insulators, which are materials 
insulating in the bulk, but which possess dissipationless 
conducting states at their edge (two-dimensional topo- 
logical insulators) or surface (three-dimensional topo- 
logical insulators). 1 ' 2 Since the introduction of the con- 
cept of two-dimensional topological insulators — often re- 
ferred to as quantum spin Hall (QSH) insulators — and 
their first prediction in graphene, 3 ' 4 several other systems 
have been proposed theoretically to exhibit QSH states, 
such as inverted HgTe/CdTe quantum-well structures, 5 
GaAs under shear strain, 6 two-dimensional bismuth, 7 or 
inverted InAs/GaSb/AlSb Type-II semiconductor quan- 
tum wells. 8 Experimentally, the QSH state has first been 
observed in inverted HgTe quantum wells, 9-12 where one 
can tune the band structure by fabricating quantum wells 
with different thicknesses. 8 Similarly to the quantum 
Hall (QH) state, which can be characterized by Chern 
numbers, 13,14 the QSH state can also be described by a 
topological invariant, in this case the invariant. 3 ' 15 
This invariant describes whether one deals with a trivial 
insulator, that is, an insulator without edge states pro- 
tected by time-reversal symmetry, or a QSH insulator. 
One of the most prominent features of QSH insulators is 
the existence of dissipationless helical edge states, that is, 
edge states whose spin orientation is determined by the 
direction of the electron momentum and are protected 
from backscattering. 16 ' 17 Thus, at a given edge, one can 
find a pair of countcrpropagating, spin-polarized edge 
states, a fact whose experimental verification has only 



very recently been reported. 12 Since those counterprop- 
agating, spin-polarized edge states are robust against 
time-reversal invariant perturbations such as scattering 
by nonmagnetic impurities, they are promising for ap- 
plications within the field of spintronics, 18,19 the central 
theme of which is the generation and control of noncqui- 
librium electron spin in solids. 

At the center of the QSH state are relativistic correc- 
tions, which can — if strong enough — lead to band inver- 
sion, that is, a situation where the normal order of the 
conduction and valence bands is inverted. 20,21 By fabri- 
cating HgTe quantum wells with a thickness larger than 
the critical thickness d c sa 6.3 nm, such an inverted band 
structure can be created in HgTe/CdTe quantum- well 
structures. In fact, materials with band inversion have 
been studied for some time 22 and another interesting 
prediction — different from the QSH state — has been that 
the combination of two materials with mutually inverted 
band structures can lead to the formation of interface 
states which — depending on the material parameters — 
can possess a linear two-dimensional spectrum. 23,24 

Following the observation of the QSH state in HgTe- 
based quantum wells, much effort has been invested in 
the theoretical investigation of the properties of two- 
dimensional topological insulators, their edge states, and 
possible applications. Examples include the extension of 
the low-energy Hamiltonian introduced in Ref. 5 to ac- 
count for additional spin-orbit terms due to out-of-plane 
inversion breaking in HgTe quantum wells 25 as well as 
studies on how helical edge states and bulk states interact 
in two-dimensional topological insulators. 26 The effect of 
magnetic fields on transport in inverted HgTe quantum 
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wells has been treated in Refs. 27-29. Moreover, the ef- 
fect of finite sizes on the QSH edge states in HgTe quan- 
tum wells has been investigated and it has been shown 
that for small widths the edge states of the opposite sides 
in a hnitc system can overlap and produce a gap in the 
spectrum. 30 Based on this coupling of the wave functions 
from opposite edges, a spin transistor based on a con- 
striction made of HgTe has been proposed. 31 Finite-size 
effects in topological insulators have not only been stud- 
ied for HgTe, but also in three-dimensional topological 
insulators, in particular the crossover to QSH insulators 
in thin films. 32 - 34 

Our purpose is to present a systematic study of the 
effect a perpendicular magnetic field has on the energy 
spectrum and magnetic edge states of HgTe/CdTe quan- 
tum wells (as described by the Hamiltonian introduced 
in Rcf. 5) in the normal as well as in the inverted regime. 
In particular, we present an analytical solution for the 
magnetic edge states confined by a hard- wall potential in 
the spirit of Refs. 35,36, where the problem of spin edge 
states and magnetic spin edge states in two-dimensional 
electron gases with hard walls and spin-orbit coupling has 
been solved analytically. Complementary to this proce- 
dure, we also make use of a numerical scheme based on 
the method of finite differences. Furthermore, the mag- 
netic properties of HgTe quantum wells are investigated 
within this model, again for both, the normal and in- 
verted regimes. 

The manuscript is organized as follows: Section II gives 
a short overview of the effective model used to describe 
the HgTe quantum well. In Sec. Ill, following the presen- 
tation of two methods to calculate the energy spectrum 
and eigenstatcs, an analytical and a finite-differences 
method, the evolution of QSH and QH states with in- 
creasing magnetic fields is discussed. The second part 
of the manuscript, Sec. IV, is devoted to the discussion 
of the magnetic properties of this system. Finally, the 
manuscript is concluded by a brief summary. 



II. MODEL 

Our model is based on the two-dimensional effec- 
tive Hamiltonian of HgTe/CdTe quantum wells derived 
from the Kane model by Bernevig et al. 5 This effec- 
tive 4x4 Hamiltonian captures the essential physics in 
HgTe/CdTe quantum wells at low energies and describes 
the spin-degenerate electron-like (E) and heavy hole-like 
(H) states \Ef}, \Hf), \E 4), and \H i) near the T point. 
The effect of a magnetic field B(r) can be included in 
this model by adding a Zeeman term 11 and promoting 
the components of the wave vector to operators, that is, 
fej — > TTi/H, where i denotes the in-plane coordinates x 
or y of the quantum well, 7Tj = pi + eA^r) the kinetic 
momentum operators, pi the momentum operators, A(r) 
the magnetic vector potential, and e = |e| the elementary 
charge. 

In our model, we consider a constant magnetic field 



perpendicular to the quantum well, that is, B = Be z 
with B > (throughout this manuscript). Since hard 
walls will be added in Sees. Ill A and III B to confine the 
system in the y-direction, it is convenient to choose the 
gauge 



A(r) = -Bye x , 
for which the effective Hamiltonian reads as 
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with the system parameters A, B, C, T>, and M., the mag- 
netic length l B = y/h/e\B\ = y/H/eB, the Bohr magne- 
ton ji B , and the 4x4 unity matrix 1. For the basis order 
\E t), | -Eft)) I -EH)) l-Ef I), the remaining 4x4 matrices 
are given by 
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where a x , a y , and a z denote the Pauli matrices and a g = 
diag(<j e , g b ) contains the effective (out-of-plane) g-factors 
g e and g h of the E and H bands, respectively. 

The material parameters introduced above, A, B, C, 2?, 
and M., are expansion parameters that, like g c and g h , 
depend on the quantum-well thickness d. 5,10 Thus, the 
quantum-well thickness can be used to tune the band 
structure. Here, A describes the coupling between the 
electron-like and hole- like bands, C and V describe a stan- 
dard parabolic dispersion of all bands, whereas M. and B 
determine whether the band structure is inverted or not: 
If the thickness of the quantum well is smaller than the 
critical thickness, d c sa 6.3 nm, the band structure is nor- 
mal and Ai /B < 0, while, for a quantum- well thickness 
above d c , the band structure is inverted and M/B > 0. 

In some cases, a reduced form of Eq. (2) can be used. 
For relatively strong magnetic fields, the terms quadratic 
with the kinetic momentum in Eq. (2) are small near the 
r point and can be omitted, as can the contribution from 
the Zeeman term, that is, B = V = and g o/h = 0. 27 ' 37 



III. MAGNETIC EDGE STATES 
A. Analytical solution 

In this section, we discuss the analytical solution — 
which in many ways resembles the calculation of the 
spin edge states in two-dimensional electron gases with 
spin-orbit coupling 35 — of the model system described by 
Eq. (2) for several different geometries: (i) bulk, that is, 
an infinite system, (ii) a semi-infinite system confined to 
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y > 0, and (iii) a finite strip with the width w in y- 
direction. For all these cases, we apply periodic bound- 
ary conditions in ^-direction. The confinement can be 
described by adding the infinite hard-wall potentials 



V(y) = 



in (ii) and 



V{y) 



for y > 
oo elsewhere 



for \y\ < w/2 
oo elsewhere 



(4) 



(5) 



in (iii). 

In order to determine the solutions for cases (i)-(iii), we 
first need to find the general solution to the differential 
equation given by the free Schrodinger equation 



(6) 



where ty(x,y) is a four-component spinor. By impos- 
ing the appropriate boundary conditions along the y- 
direction on this general solution, we can obtain the solu- 
tions for each of the cases considered. Since translational 
invariance along the x-direction as well as the spin direc- 
tion are preserved by H and H + V(y)l, respectively, the 
wave vector in x-direction, k, and the spin orientation, 
s =t / i, are good quantum numbers in each of the three 
cases, which naturally suggests the ansatz 





where L is the length of the strip in x-direction and 
where, for convenience, we have introduced the trans- 
formation £ = £(y) = -y/2 (y - l%k) /l B . 

Inserting the ansatz (7) for spin-up electrons into 
Eq. (6), we obtain the following system of differential 
equations: 
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Due to the specific form of Eq. (8), its solution can be 
conveniently written in terms of the parabolic cylindrical 
functions D u (£), which satisfy the following recurrence 
relations: 38 
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With the heavy hole- like component (£) coupled to the 
electron-like component /f (£) by the raising operator and 
the opposite coupling described by the lowering operator, 
one type of solution is of the form 

M0=viDu(0 and fft(0=«aA,_i(£), (11) 

where V\ and are complex numbers, which are to be 
determined by solving the system of linear equations ob- 
tained from inserting this ansatz into Eq. (8). This sys- 
tem has non-trivial solutions for 
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By determining those non-trivial solutions, for A ^ we 
find the two (non-normalized) solutions 



xl(0 = (V^ 4 (0/Jb,4A4-i(0 

to Eq. (8) with 



(15) 
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However, there is a second set of — in general — 
independent solutions to Eq. (8) that can be obtained 
from the ansatz 

/ t (0=«iA,(-0 and fft(0 = «2A,-i(-f). (17) 

where u\ and u-i are complex numbers as before. With 
this ansatz yielding two further solutions, 



V l(0 = (V2AD 1/1 (-0/l B ,-clD vl _ 1 (-0) , (18) 
the general solution to Eq. (8) — if A ^ — is given by 



( M0 \ 

\9t(0J 



!xl(Q+0xU®+ivl(Q+8vHO, (19) 



where the coefficients a, ft, 7, and 5 are complex num- 
bers to be determined by the boundary conditions of the 
problem. 

A procedure similar to the one above can also be ap- 
plied for the spin-down electrons in Eq. (7). Then, we 
find 



\9l(0 



= axX(C)+P X 1 - (0 + 7 Vi (0 + 6 rt (0 , (20) 
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where we have introduced the vectors 

xt(0 = (ciD pi _^),V2AD ui ^)/l B )' 



and 
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and 
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_(24) 

As in the case of spin-up electrons, the coefficients a, /3, 
7, and 5 need to be fixed by boundary conditions. In 
the following, we will use the general solutions given by 
Eqs. (19) and (20) to determine the energy spectrum and 
wave functions for several different geometries. 
(i) Bulk. 

If there is no confining potential V(y), that is, if 
we consider an infinite system, where Eq. (8) holds for 
any £ £ R, we only have to require the wave func- 
tion to be normalizable and accordingly we impose the 
boundary conditions lim /t(£) = lim g-\(£,) = 0. 

>±oo £— >±oo 

These requirements can only be satisfied if v is a non- 
negative integer n in Eq. (11). In this case, D n (£) = 
2~ n / 2 Q-6 / 4 7J TI (£/v / 2) can be expressed by Hcrmite poly- 
nomials i? ra (£), 38 an d both Eqs. (11) and (17) lead to the 
same solution. If n > 1, the ansatz from Eq. (11) leads 
to an eigenvalue problem for E from which the following 
Landau levels for spin-up electrons can be determined: 
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For n = 0, on the other hand, Eqs. (11) and (17) reduce 
to the ansatz /f(£) = ^iA)(0 an d <7t(£) = an d we 
obtain the Landau level 

E^(0) = C + M-^^- + ^ B B. (26) 
By requiring lim = lim <?^(£) = 0, the Lan- 

£— ^±oo £— f±oo 

dau levels for spin-down electrons can be calculated sim- 
ilarly as 
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With Eqs. (25)- (28), we have recovered the Landau levels 
found in Rcf. 11. The corresponding cigenstates are given 
in the Appendix A. 

In writing down Eqs. (25)- (28), we have adopted the 
convention that B > 0, that is, the magnetic field points 
in the z-direction. The formulas of the Landau lev- 
els for B < can be obtained from Eqs. (25)-(28) via 
the relations E^O.B) = E~ s (0,-B) and E±(n, B) = 
E± s (n, —B) [note that the magnetic length in Eqs. (25)- 

(28) is given by l B = y/h/e\B\]. 
(ii) Semi- infinite system. 

In the presence of the confining potential given by 
Eq. (4), the wave function is required to vanish at the 
boundary y = as well as at y — > oo. Thus, we invoke 
the boundary conditions lim /-wi (£) = lim <7 T /i(£) = 

and /f/4.(£o) = 9f-/±(€o) — for spin- up as well as spin- 
down electrons, where £o = —\2l B k. The condition for 
£ —> oo can only be satisfied for 7 = 5 = and 7 = 5 = 0, 
respectively. Then, each remaining pair of coefficients, a 
and /J as well as a and /?, from Eqs. (19) and (20) has to 
be calculated from the condition at y = 0, that is, at £o- 
The resulting linear systems of equations have non-trivial 
solutions if 



J/l 



- c y i D vfn _ 1 (Z )D vVl (to) = 0. 



(29) 



This transcendental equation enables us to calculate the 
electron dispersion for spin-up [s =f in Eq. (29)] as well 
as for spin-down electrons [s =1 in Eq. (29)]. The cor- 
responding eigenstates can be determined by explicitly 
calculating the coefficients a, /3 and 5,/3, respectively. 
(Hi) Finite-strip geometry. 

In the finite-strip geometry described by Eq. (5), the 
wave function has to vanish at the potential bound- 
aries, that is, Eqs. (19) and (20) have to vanish at 
£1/2 = V2 (tw/2 — l 2 k) jl B . The corresponding linear 
systems of equations defined by this condition have non- 
trivial solutions if 



xl A (£i) x^(a) v¥*(ti) 

xl A (6) xl A (6) rftHb) vt'Hb) 

(30) 

for spin-up (s = J [) and spin-down (s =^) electrons, re- 
spectively. Similarly to (ii), the transcendental Eq. (30) 
represents exact expressions from which the dispersion 
of the electrons can be calculated. The corresponding 
eigenstates can be determined by explicitly calculating 
the coefficients a, (3, 7, and 5 for spin-up electrons and 
a, /3, 7, and 5 for spin-down electrons, respectively. 

Having derived transcendental equations from which 
the electronic dispersion (and indirectly the eigenstates) 
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can be determined for semi-infinite as well as finite-strip 
systems, we will also introduce an alternative method to 
calculate the spectrum and eigenstates of a finite strip. 



B. Numerical finite-difference solution 

In addition to solving the exact expression (30), we 
calculate the eigenspcctrum and eigenstates also by us- 
ing a finite-difference scheme to express Eq. (2). 39 We 
discretize Eq. (2) for B = and account for the mag- 
netic field by introducing the Peierls' phase 40 to describe 
the vector potential given by Eq. (I) and an additional 
on-site term to describe the Zeeman term. If only nearest 
neighbors are considered and there is no magnetic field, 
this procedure leads to the Hamiltonian introduced in 



Ref. 10. 

For reasons of improving the convergence of our calcu- 
lation, we go beyond the nearest-neighbor approximation 
and include the next-nearest neighbors. Due to transla- 
tional invariance along the x-direction, the ^-coordinate 
can be Fourier transformed to the reciprocal space and 
wc obtain the Hamiltonian 



FD 
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(31) 



a/3 



where k is the momentum along the x-direction, n and 
n' G Z are discrete y-coordinates, a and /3 denote the 
basis states \Ef), \H"[), \E I), \H I), and c\ na (c kna ) 
is the creation (annihilation) operator of those states. 
Furthermore, we have introduced the matrix 
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3 6 
and a denotes the distance between two lattice points in 
y-direction. However, in the finite-strip geometry con- 
sidered here, the matrix given by Eq. (32) has to be 
modified at the edges along the y-direction, where only 
nearest neighbors can be used for the approximation of 
the derivatives with respect to y. Following these mod- 
ifications, the eigenspcctrum and the eigenstates of the 
system in a finite-strip geometry can be determined nu- 
merically. 



Comparison between the analytical and 
numerical solutions 



We compare the results obtained by the analytical pro- 
cedures described in Sec. Ill A with those of the finite- 
difference method introduced in Sec. IIIB. For illustra- 
tion, Fig. 1 shows the energy spectra of a semi-infinite 



system [Fig. I (a)] and a finite strip of width w = 200 
nm [Fig. I (b)]. Here, wc have chosen the magnetic 
field B = 10 T and the parameters A = 364.5 meV 
nm, B = -686.0 meV nm 2 , C = 0, V = -512.0 meV 
nm 2 , M = —10.0 meV, and g c = g h = 0, which (apart 
from the vanishing g-factors) correspond to the thickness 
of d = 7.0 nm. 2,10 Whereas the energy spectrum of a 
semi-infinite system is calculated using the transcenden- 
tal Eq. (29), both procedures described above, solving 
the transcendental Eq. (30) or diagonalizing the finite- 
difference Hamiltonian (31), can be used to calculate the 
eigenspectrum of the Hamiltonian (2) in a finite-strip ge- 
ometry. The finite-difference calculations for Fig. 1 (b) 
have been conducted for 201 lattice sites along the y- 
direction, for which we get a relative error of 10 -6 -10~ 5 . 
Figure 1 (b) also clearly illustrates the nearly perfect 
agreement between the analytical and numerical solu- 
tions. As can be expected if the magnetic length l B is 
small compared to the width of the sample w, the energy 
spectra near the edge as well as the energy spectra in the 
bulk are almost identical for the semi-infinite and finite 
systems as shown in Figs. 1 (a) and 1 (b). The bulk Lan- 
dau levels are perfectly characterized by Eqs. (25)- (28). 
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y k N y k [nm] 

FIG. 1: (Color online) Calculated energy spectra of (a) a 
semi-infinite system and (b) a finite strip of width w = 200 
nm for B = 10 T, A = 364.5 meV nm, B = -686.0 meV 
nm 2 , C = 0, V = -512.0 meV nm 2 , M = -10.0 meV, and 
9c — <7h — 0. Here, the energy spectra are plotted versus 
yk = l%k. The solid and dashed lines represent s =\ and 
s =4- states, respectively, which have been calculated using 
the analytical methods from Sec. Ill A [case (ii) for Fig. (a) 
and case (iii) for Fig. (b)]. Results obtained by the finite- 
difference method from Sec. IIIB are represented by circles 
(spin up) and diamonds (spin down) in Fig. (b). 
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FIG. 2: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = \^(x,y)\ 2 , of selected 
states for d = 5.5 nm, w = 200 nm, and B = T, where solid 
and dashed lines represent s =t and a =4- states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the x-direction is given by 
v k = [dE(k)/dk]/h. 
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FIG. 3: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = \"^(x,y)\ 2 , of selected 
states for d — 5.5 nm, w = 200 nm, and B = 0.1 T, where 
solid and dashed lines represent s =f and s =4- states, respec- 
tively. Here, the states shown in Figs, (b) and (c) are marked 
in the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the a-direction is given by 
v k = [dE(k)/dk\/h. 
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FIG. 4: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = \^(x,y)\ 2 , of selected 
states for d = 5.5 nm, w = 200 nm, and B = 1 T, where solid 
and dashed lines represent s and s =4- states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the i-direction is given by 
v k = [dE(k)/dk]/h. 



also Sec. IIIC). 



D. Results 

In this section, we investigate the magnetic field de- 
pendence of the energy spectrum and its correspond- 
ing eigenstates in a finite-strip geometry with the width 
w = 200 nm. The graphs shown in this section have 
been calculated using the finite-difference scheme from 
Sec. IIIB with 201 lattice sites along the y-direction (see 



1. Ordinary insulator regime 

First, we examine the quantum- well spectrum in the 
ordinary insulator regime, that is, for a thickness d < 
d c , where the band structure is normal and there are no 
QSH states (at zero magnetic field). Figures 2-5 show the 
energy spectrum and (selected) eigenstates at different 
magnetic fields for the material parameters A = 387 meV 
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FIG. 5: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = \^(x,y)\ 2 , of selected 
states for d = 5.5 nm, w = 200 nm, and B = 10 T, where solid 
and dashed lines represent s =f and s =4- states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the ^-direction is given by 
v k = [dE(k)/dk]/h. 



FIG. 6: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = | \& (a;, j/) | 2 , of selected 
states for d = 7.0 nm, w = 200 nm, and B = 0T, where solid 
and dashed lines represent s =j" and s —\. states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the a-direction is given by 
v k = [dE(k)/dk]/h. 
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nm, B = -480.0 meV nm 2 , C = 0, V = -306.0 mcV nm 2 , 
and M. = 9.0 meV, which correspond to a quantum-well 
thickness of d = 5.5 nm. 2 As illustrated by Fig. 2 (a), 
which shows the spectrum for B = 0, only bulk states, 
but no edge states can be found [see Figs. 2 (b) and (c)], 
a situation which changes little if small magnetic fields 
are applied (see Fig. 3). Only if the magnetic field is 
increased further, do Landau levels [given by Eqs. (25)- 
(28)] and corresponding QH edge states begin to form as 
can be seen in Figs. 4 and 5. Comparing Figs. 4 and 5, 
one can also discern that with increasing magnetic field 
the QH edge states become more localized. 
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2. QSH regime 

In Figs. 6-9, by contrast, the energy spectrum and (se- 
lected) eigenstates of a strip with the width w = 200 
nm are presented for the material parameters A = 364.5 
meV nm, B = -686.0 mcV nm 2 , C = 0, V = -512.0 mcV 
nm 2 , M = -10.0 meV, g e = 22.7, and g h = -1.21, cor- 
responding to a quantum-well thickness d = 7.0 nm, ' 10 
that is, for parameters in the QSH regime (at B = 0), 
and several strengths of the perpendicular magnetic field. 
The spectra and states in Figs. 6-9 illustrate the evolu- 
tion of QSH and QH states in HgTe. 

Figure 6 (a) shows the spectrum at zero magnetic field. 
At this magnetic field, one can observe the QSH state in- 
side the bulk gap, that is, two degenerate pairs of coun- 
terpropagating, spin-polarized edge states, one pair at 
each edge [see Figs. 6 (b) and (c)]. As found in Ref. 30, 
at k = the wave functions of QSH edge states with the 
same spin, but at opposite edges overlap thereby open- 
ing up a gap [see the inset in Fig. 6 (a)]. By increasing 



FIG. 7: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = |>l/(:r,2/)| 2 , of selected 
states for d — 7.0 nm, w = 200 nm, and B = 0.1 T, where 
solid and dashed lines represent s =t and s =4- states, respec- 
tively. Here, the states shown in Figs, (b) and (c) are marked 
in the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the i-direction is given by 

v k = [dE(k)/dk]/h. 

the width of the strip, the overlap of the edge-state wave 
functions with the same spin is diminished and one can 
remove this finite-size effect. 

For small magnetic fields (Fig. 7), apart from the split- 
ting of spin-up and down states, the situation is at first 
glance quite comparable to the one in Fig. 6. Most im- 
portantly, one can still find pairs of counterpropagat- 
ing, spin-polarized states in the vicinity of each neutral- 
ity point [for example, the states shown in Figs. 7 (b) 
and (c)], that is, the crossovers between the lowest (hole- 
like) conduction band and uppermost (electron-like) va- 
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FIG. 8: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = | \& (2c, j/) | 2 , of selected 
states for d = 7.0 nm, w = 200 nm, and B = 1 T, where solid 
and dashed lines represent s ="|" and s =J, states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the x-direction is given by 
v k = [dE(k)/dk]/h. 



FIG. 9: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = ^(a;, y) | 2 , of selected 
states for d = 7.0 nm, w = 200 nm, and B — 10 T, where solid 
and dashed lines represent s =j" and s =4, states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the i-direction is given by 
v k = [dE(k)/dk]/h. 



lence band [marked by dots in Fig. 7 (a)]. However, 
we stress that these counterpropagating, spin-polarized 
states which can be found (at a given edge) if the Fermi 
level is close to the neutrality points, are not connected 
with each other by time-reversal symmetry and are there- 
fore not topologically protected (for example, against 
spin-orbit coupling). 

Going to B = 1 T (Fig. 8), we can still find coun- 
terpropagating, spin-polarized states near and at the 
crossovers between the lowest (hole-like) conduction and 
uppermost (electron-like) valence bands, which (in the 
bulk) have evolved into the E^(Q) and E^(0) Landau 
levels. As the center of the orbital motion is given by 
one can sec that those states are now no longer 
as localized as before at the edges [see Figs. 8 (b) and (c)]. 
Meanwhile, the bulk states from Fig. 6 have also evolved 
into Landau levels given by Eqs. (25) and (27) with lo- 
calized QH edge as well as bulk states. From Fig. 8, one 
can also discern another feature of the energy spectrum 
and eigenstatcs that develops with an increasing mag- 
netic field, namely the appearance of 'bumps' [see the 
spin- up valence bands in Fig. 8 (a)]. If the Fermi level 
crosses those 'bumps', one finds states which are localized 
near the same edge and carry the same spin, but counter- 
propagate. This has also been observed in Ref. 29, where 
those states gave rise to exotic plateaus in the longitudi- 
nal and Hall resistances. As can be seen in Figs. 4 and 5 
(as well as later in Figs. 9, 14, and 15), this behavior can 
also be found for other quantum well parameters. 

The situation described so far changes for high mag- 
netic fields (Fig. 9), when the electron-like band de- 
scribed by -E^(O) (in the bulk) is above the hole- like 
-E^(O) band. Then, there is no longer any crossover be- 
tween the dispersions of electron- and hole- like bands and 



one consequently cannot find counterpropagating, spin- 
polarized states anymore, just QH edge states propagat- 
ing in the same direction [for example, the states shown 
in Figs. 9 (b) and (c)]. 

As has been known for a long time, the uppermost 
(electron-like) valence and the lowest (hole-like) conduc- 
tion Landau levels cross at a finite magnetic field B c in 
inverted HgTe/CdTe quantum wells. 41-43 The transition 
between the two situations, the one where counterprop- 
agating, spin-polarized states exist and the one where 
they do not, happens exactly at this crossover point: As 
long as the hole-like band is above the electron-like band, 
that is, as long as the band structure remains inverted, 
one can find counterpropagating, spin-polarized states in 
addition to the QH states. Otherwise, there are only QH 
states. 

This crossover point can be easily calculated from the 
Landau levels via the condition E^(0) = E^(0), from 
which we get 

M 

Bc = 2nB/* -(g c +g b )» B /4 (35) 

for the magnetic field at which the transition happens 
(valid only for B c > 0). Here, $0 = 27rft/e denotes the 
magnetic flux quantum. The validity of the result given 
by Eq. (35) is also illustrated by Fig. 10, which shows 
the magnetic field dependence of the energies of the finite 
strip with width w — 200 nm at k — and of the bulk 
Landau levels for the same band parameters as above. 
As can be expected, the energies at k = are given 
by the Landau levels (25)- (28) at high magnetic fields. 
Most importantly, the crossover between the electron- 
like -E^(O) and the hole- like £^(0) bands happens in the 
region, where the independence of the energy levels at 
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B [T] 

FIG. 10: (Color online) Magnetic field dependence of the 
states at k — in a finite strip of width w = 200 nm com- 
pared to the bulk Landau levels given by Eqs. (25)-(28). The 
thinner solid and dashed lines represent bulk Landau levels 
for s =t and s =4-, respectively. The levels of the finite-strip 
geometry are displayed by thick lines. All levels displayed 
here have been calculated for band parameters corresponding 
to d = 7.0 nm. 
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FIG. 12: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = | \& (a;, j/) | 2 , of selected 
states for d = 6.3 nm, w = 200 nm, and B = 0T, where solid 
and dashed lines represent s =j" and s =4- states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the a-direction is given by 
v k = [dE(k)/dk]/h. 
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FIG. 11: (Color online) Magnetic field dependence of the 
states at k = in finite strips with the widths (a) w — 25 nm, 
(b) w — 50 nm, (c) w = 75 nm, and (d) w = 100 nm com- 
pared to the bulk Landau levels given by Eqs. (25)-(28). The 
thinner solid and dashed lines represent bulk Landau levels 
for s =t and s =4-, respectively. The levels of the finite-strip 
geometry are displayed by thick lines. All levels displayed 
here have been calculated for band parameters corresponding 
to d = 7.0 nm. 



in Sec. II) for HgTe has been used, and Ref. 29. Only for 
high magnetic fields, the band structure becomes nor- 
mal and one enters the ordinary insulator regime, in 
which no countcrpropagating, spin-polarized states can 
be found (see also Ref. 29). We remark that the de- 
scription presented in this section also bears out if other 
widths w > 100 nm of the finite strip are investigated. 
For larger widths, the formation of Landau levels sets in 
already at lower magnetic fields, whereas higher fields are 
needed to observe Landau levels in more narrow strips. 
If very small samples (w < 50 nm) are investigated, 
however, we find that there is no crossover between the 
electron-like £^(0) and the hole-like £^(0) bands, as il- 
lustrated by Fig. 11, which shows a comparison between 
the bulk Landau levels and the states calculated at k = 
for band parameters corresponding to d = 7.0 nm and 
several small widths w. Only if w > 50 nm, the gap due 
to the finite size of the sample at B = is reduced far 
enough and one can observe a crossover of the (0) and 
£^(0) bands at B = B c which is then give by Eq. (35). 



3. Critical regime 



k = is already described extremely well by those Lan- 
dau levels and from Eq. (35) we find B c pa 7.4 T, con- 
sistent with the numerical result that can be extracted 
from Fig. 10. Furthermore, one can see how the £^(0) 
band is below the £7*(0) band for B < B c , and how the 
situation is reversed for B > B c . 

Therefore, we find that if the magnetic field is not too 
high, the countcrpropagating, spin-polarized states per- 
sist at finite magnetic fields, consistent with the conclu- 
sions in Refs. 27,28, where the reduced model (mentioned 



Finally, for the purpose of comparison to the discussion 
above, Figs. 12-15 show the energy spectrum and (se- 
lected) eigenstates at different magnetic fields for a strip 
with the width w = 200 nm and the material parame- 
ters A = 373.5 meV nm, B = -857.0 meV nm 2 , C = 0, 
V = -682.0 meV nm 2 , M = -0.035 meV, g c = 18.5, 
and g h = 2.4, which correspond to the critical regime at 
a quantum-well thickness of d = d c — 6.3 nm. 2,11 For 
B = 0, instead of edge states, we find states whose prob- 
ability densities are spread over the entire width of the 
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FIG. 13: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = | \& (2c, j/) | 2 , of selected 
states for d — 6.3 nm, w = 200 nm, and B = 0.1 T, where 
solid and dashed lines represent s =f and s =4- states, respec- 
tively. Here, the states shown in Figs, (b) and (c) are marked 
in the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the s-direction is given by 
v k = [dE(k)/dk]/h. 
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FIG. 15: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = | \& (2c, j/) | 2 , of selected 
states for d = 6.3 nm, w = 200 nm, and B — 10 T, where solid 
and dashed lines represent s =j" and s =4, states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the i-direction is given by 
v k = [dE(k)/dk]/h. 



IV. MAGNETIC OSCILLATIONS 
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FIG. 14: (Color online) (a) Calculated energy spectrum and 
(b), (c) probability densities, p(x,y) = \^(x,y)\ 2 , of selected 
states for d = 6.3 nm, w = 200 nm, and B = 1 T, where solid 
and dashed lines represent s =f and s =\. states, respectively. 
Here, the states shown in Figs, (b) and (c) are marked in 
the energy spectrum, Fig. (a), by dots. The velocity with 
which the states propagate along the ^-direction is given by 
v k = [dE{k)/dk]/h. 



strip with a slight preponderance near one of the edges 
[see Figs. 12 (b) and (c)]. with increasing magnetic field 
the states become more localized (see Figs. 13 and 14) 
and, finally, one can find QH edge states (see Fig. 15). 



A. General formalism 

In this section, we discuss the magnetization and mag- 
netic oscillations in HgTe quantum wells. Our starting 
point is the grand potential 



n{T,n,B) 



dep(e) ln{l + exp[-/3(e- 



(36) 

where /3 = l/(k B T) and T denotes the temperature, k B 
the Boltzmann constant, /1 the chemical potential, p(e) 
the density of states per unit area, and S is the surface 
area. 

We make the electron-hole transformation and divide 
the spectrum in the electron and hole contributions, 
p e (e) = p(e)Q(e - E n ) and p h (e) = p(e)Q(E n - e), where 
E n = E n (B) denotes the neutrality point. Then, we can 
rewrite f2 (T, fi, B) as 



n (t, n, b) =o e (r, /i, b) + n h (T, n, B) 

+ S de ph(e) (e - n) , 



(37) 



where 



a {T,n,B) = ~ J dep c (e)ln{l + exp[-/3( £ - A1 )]} 

(38) 



and 

n h (T,fi,B) 



s 



dep b (e) ln{l + cxp [0 (e-ft)]} 

(39) 
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denote the grand potentials of electrons and holes, re- 
spectively. The total particle number in the system is 
given by N tot = — [dfl (T, ll, B) /dfj]. However, it is more 
convenient to distinguish between electrons and holes and 
to work with the carrier imbalance N = N — N h (with 
iV e / h denoting the number of electrons and holes, respec- 
tively). Following Ref. 44, we redefine the grand potential 
and use 

Q' (T, (i, B) = fi (T, fi, B) + Sfi J de p h (e) 

= a (t, fx, b) + n h (T, n, B) + n (B), 



where 



ft {B) = S J dep h {e)t 



(40) 



(41) 



is the ground-state/ vacuum energy. The carrier imbal- 
ance is then given by N = — [dCl' (T, fi, B) /d/j]. 

The magnetization (as a function of the chemical po- 
tential, the temperature, and the magnetic field) can be 
extracted from CI' (T, \i, B) via 

IdW (T,ti,B) 
S~ 



M tot (T,li,B) 



dB (42) 
= M (B)+M (T,fi,B), 

where we have split the magnetization in the vacuum 
part 

1 dQ (B) 
~S dB 

and the non- vacuum part 



M Q {B) 



(43) 



M(T,fx,B) = -^ 



dn e (T,fi,B) + dfl h (T,fi,B) 



dB 



dB 



(44) 

At zero temperature, the magnetization of an undoped 
system is given by Mq(B), whereas at finite tem- 
peratures or in doped systems the additional contri- 
bution M (T, li, B) arises. The magnetization as a 
function of the carrier imbalance density rid = N / S 
(rid > : n-doped, < : p-doped) is given by 
M [T, [i (T, rid, B) , B] , where the chemical potential is de- 
termined by 



1 

S 



dVl' (T, n, B) 



dfx 



(45) 



p, =l i(T,n d ,B) 



Finally, we remark that the magnetic susceptibility 
Xtot (T,fi,B) = xo(B) + x {T,fi,B) can also be split in 
the vacuum part 

dM Q (B) _]_d 2 n (B) 



Xo(B) = 



dB 



and the non-vacuum part 



dB 2 



dM h (T,n,B) 



(46) 



dB 

d 2 n c (T,Li,B) 

dB 2 



dB 

d 2 n h (T jt i,B) 

dB 2 



(47) 



For the (bulk) Landau levels (and typical parameters 
of HgTc quantum wells), the different contributions to 
the grand potential read as 

fl Q (T,n,B) = 



SB 
Wo 



In 



l+e-^^W-f)] e [^(o)-^(0)] 



+ In [l + e-^( 0) -")] 9 [EHO) - ^(0)] 



J2 E ln[l + e^( £ + (n) -^ 

n=l s=t>4 



(48) 



Vl h (T,n,B) = 
SB ' 



/3$ c 



In 



1 + e 



9[£ T (0)-£^(0)] 



+ In [l + e ^( st (°)-^)] 9 [E l (0) - £^(0)] 



and 



n=l s=t,i 



n (B) = n dls (B) + n (B) 



(49) 



(50) 



where the energies are given by Eqs. (25)-(28) and 4>o = 
2irh/e is the magnetic flux quantum. In Eq. (50), we have 
split the ground-state potential into a contribution from 
the uppermost valence band [which may not be contin- 
uously differentiable if there is a crossover between the 
hole-like E^(0) and the electron-like E^(0) bands like at 
the transition point in Fig. 10], 



fl dis (B) =£+(0)9 [E^(0) - E l (0)] 

+ E t (0)Q[E i (0)-E*(0)] , 
and a contribution from the remaining valence bands, 



(51) 



= E ]T Ei(n). 



(52) 



n—1 • 



Since the energies in Eq. (52) are not bounded from 
below (for typical parameters of HgTc quantum wells), 
the sum is divergent; following Refs. 45-47, we introduce 
a smooth cutoff function which results in a smooth rio(B) 
(we refer to the Appendix B for more details). If there 
is no crossover between the electron-like Er (0) band and 
the hole-like E^(Q) band, that is, if one deals with an 
ordinary insulator, then the total ground-state magneti- 
zation Mo(B) is continuous. Due to fl dis (B), which is not 
continuously differentiable if the £^(0) and £^(0) bands 
cross (see Fig. 10), the ground-state magnetization is not 
continuous at the crossover point in this case. For bulk 
Landau levels, we find the jumps 

AM = lim [M (B C + 5B) - M (B C - SB)} = - — 

(53) 
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FIG. 16: (Color online) The non-vacuum magnetization 
M (T, n, B) (corresponding to a quantum-well thickness of 
d = 7.0 nm) plotted versus l/B for a fixed chemical potential, 
(i = 20 meV, and different temperatures (T = 1, 10, 100 K). 
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FIG. 17: (Color online) The non-vacuum susceptibility 
X (T, (j,, B) (corresponding to a quantum-well thickness of 
d = 7.0 nm) plotted versus l/B for a fixed chemical potential, 
fi — 20 meV, and different temperatures (T = 1, 10, 100 K). 



at the crossover point B c , where there is a transition 
from the inverted [E*(0) < E^(0)] to the normal regime 
[£^(0) < E^(0)}. 

However, at finite temperatures or doping, the total 
magnetization is given by the sum of the ground-state 
magnetization Mq(B) and the contribution from the elec- 
trons and holes, M (T, /i, B). Analyzing this contribu- 
tion for the case of a transition from the inverted to the 
normal band structure, one finds that M (T, \i, B) van- 
ishes for zero temperature and zero doping, but other- 
wise always contains a discontinuity at B c which exactly 
cancels the discontinuity of the intrinsic magnetization. 
Thus, the total magnetization is a continuous function. 
If there is no transition between the normal and inverted 
band structures, the non- vacuum contribution and there- 
fore the total magnetization are also continuous. For a 
given quantum- well thickness d, the vacuum contribution 
Mq{B) constitutes the same background for every set of 
thermodynamic variables (/x, T) or (rid, T) of the sys- 
tem. Thus, the quantity of interest which allows one to 
compare different doping levels, chemical potentials or 
temperatures of the system is the non- vacuum contribu- 
tion M (T,fi,B). 

Equations (43)-(50) allow us to calculate the (bulk) 
magnetization and susceptibility of HgTc quantum wells, 
the results of which are discussed in the following section. 



B. Results 

In this section, we apply the formalism introduced 
above to calculate the bulk magnetization of HgTc for the 
parameter set corresponding to a quantum- well thickness 
of d = 7.0 nm (nominally the QSH regime; see above), 
that is, a situation where there is a crossover between the 
£f(0) and £^(0) bands. Figures 16 and 17 show the mag- 
netic field dependence of the non- vacuum contributions, 



that is, the contribution arising from electrons and holes, 
to the magnetization and the susceptibility for a fixed 
chemical potential, several different temperatures, and 
magnetic fields well below the crossover point B c pa 7.4 
T (compare to Sec. HID). As different Landau levels 
cross the Fermi level with increasing magnetic field, one 
can observe the de Haas-van Alphen oscillations in the 
magnetization as well as in the susceptibility whose am- 
plitude decreases with increasing temperature. For high 
magnetic fields (see the inset in Fig. 17), the spacing 
between the energies of spin-up and spin-down Landau 
levels (with the same quantum number n) is large enough 
compared to thermal broadening to observe spin-resolved 
peaks in the susceptibility. Fitting the oscillations of the 
magnetization to a periodic function, we find that the pe- 
riodicity of those oscillations is given by A(l/B) ps 1.43 
1/T [see also the Appendix C, where Eq. (CI 7) yields a 
period of A (l/B) ps 1.35 1/T for the main contribution 
to the oscillations in the reduced model] . 

Next, we consider a fixed carrier density > 0. The 
corresponding chemical potential as a function of the 
magnetic field is calculated via Eq. (45) and is displayed 
in Fig. 18 for the density rid = 10 16 1/m 2 and different 
temperatures. With varying magnetic field, the Fermi en- 
ergy (i (0, rid, B) shows oscillations consisting of a pair of 
spin-resolved peaks, where each of those oscillations cor- 
responds to a crossing of a Landau level with the Fermi 
level. Higher temperatures result in a smoothening of the 
oscillations and a diminution of their amplitudes. More- 
over, thermal broadening leads to a removal of the spin- 
resolution at small magnetic fields. 

Figures 19 and 20 show the chemical potential and the 
combined contribution M dia (B) + M (T, n, B) to magneti- 
zation as functions of the magnetic field for T = 10 K and 
different carrier densities rid- {Here, we have added the 
discontinuous contribution from the ground-state mag- 
netization, M diB (B) = -(l/S)[dn dlB (B)/dB], to the non- 
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FIG. 18: (Color online) Magnetic field dependence of the 
chemical potential \x (T, rid,B) (corresponding to a quantum- 
well thickness of d — 7.0 nm) for rid = 10 16 1/m 2 and different 
temperatures (T = 0, 10, 100 K). 
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FIG. 20: (Color online) Magnetic field dependence of the con- 
tribution M di s(B) + M (T, fi, B) (corresponding to a quantum- 
well thickness of d = 7.0 nm) for T — 10 K and different 
densities (n d = 10 14 , 10 15 , 10 16 1/m 2 ). 
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FIG. 19: (Color online) Magnetic field dependence of the 
chemical potential \x (T, rid,B) (corresponding to a quantum- 
well thickness of d = 7.0 nm) for T = 10 K and T = 100 K 
and different densities (n d = 10 14 , 10 15 , 10 16 1/m 2 ). 
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FIG. 21: (Color online) The non-vacuum magnetization 
M (T, n, B) (corresponding to a quantum-well thickness of 
d = 5.5 nm) plotted versus 1/B for a fixed chemical potential, 
/i = 20 meV, and different temperatures (T = 1, 10, 100 K). 



vacuum magnetization in order that the discontinuity at 
B = B c be canceled.} As above, one can see the de 
Haas-van Alphen oscillations in the magnetization (see 
Fig. 20), which — for fixed carrier densities — follow the 
oscillations in the chemical potential (see Fig. 19). At 
low densities, on the other hand, only the lowest conduc- 
tion Landau level is occupied and the chemical potential 
roughly follows this level and there are consequently no 
oscillations. 

For the sake of comparison to the situation in the in- 
verted regime discussed so far, Figs. 21 and 22 show the 
magnetic field dependence of the non- vacuum contribu- 
tions to the magnetization and the susceptibility in the 
normal regime (corresponding to the parameters for a 
quantum- well thickness of d = 5.5 nm as in Sec. Ill D) 
for a fixed chemical potential and several different tem- 
peratures. As in Figs. 16 and 17, one can observe the de 
Haas- van Alphen oscillations. No discernible features are 



seen when comparing the inverted and normal regimes in 
the bulk. 

In limiting cases, compact analytical formulas to de- 
scribe some of the main features of the magnetization 
and the susceptibility shown above can be given for the 
reduced model and are presented in the Appendix C. 



V. CONCLUSIONS 

We have derived analytical formulas to calculate the 
energy spectra of HgTe quantum wells in infinite, semi- 
infinite, and finite-strip systems in the presence of per- 
pendicular magnetic fields and hard walls. Complemen- 
tary to the analytical formulas, we have also used a 
finite-difference scheme to investigate the magnetic field 
dependence of the energy spectra and their respective 



14 




-8n i i i i i i i 
0.25 5 10 15 20 

1/B [1/T] 



densities follow the oscillations in the chemical potential. 
Corresponding to those oscillations of the magnetization, 
on can also observe oscillations in the magnetic suscep- 
tibility. With increasing temperature, the amplitude of 
these oscillations decreases. Furthermore, we found that, 
if the band structure is inverted, the ground-state magne- 
tization (and consequently also the ground-state suscep- 
tibility) is discontinuous at the crossover point between 
the uppermost valence and lowest conduction Landau 
levels. At finite temperatures and/or doping, however, 
this discontinuity is canceled by the contribution from 
electrons and holes and the total magnetization and sus- 
ceptibility are continuous. 
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FIG. 22: (Color online) The non- vacuum susceptibility Tmg WQrk wafj supported by thc DFG SFB 689 and 

X\T, n, B) (corresponding to a quantum-well thickness of Gf{_K 1570 
d = 5.5 nm) plotted versus 1/B for a fixed chemical potential, 
fi = 20 meV, and different temperatures (T = 1, 10, 100 K). 



cigenstatcs in a finite-strip geometry for parameters cor- 
responding to the normal (d < d c ), inverted (d > d c ), 
and critical regimes (d ~ d c ). In the inverted regime 
(d > d c ), we found that for magnetic fields below the 
crossover point between the uppermost (electron-like) va- 
lence and lowest (hole-like) conduction Landau levels, 
one can still observe counterpropagating, spin-polarized 
states at finite magnetic fields, although these states are 
no longer protected by time-reversal symmetry Above 
the crossover point, the band structure becomes normal 
and one can no longer find those states. This situation 
is similar for parameters corresponding to the normal 
regime (d < d c ), where one cannot find counterpropagat- 
ing, spin-polarized states even for zero or weak magnetic 
fields. Finally, we have studied the bulk magnetization 
and susceptibility in HgTe quantum wells and have in- 
vestigated their dependence on the magnetic field, chem- 
ical potential, and carrier density. In the case of fixed 
chemical potentials as well as in the case of fixed den- 
sities, the magnetization (for both, the normal as well 
as the inverted regime) exhibits characteristic de Haas- 
van Alphen oscillations, which in the case of fixed carrier 



Appendix A: Landau levels 

In the absence of any confining potential, we require 
the wave functions given by Eqs. (19) and (20) to vanish 
for £ — > ±oo, which can only be satisfied if thc indices 
of the parabolic cylindrical functions are non-negative 
integers n. As above, we first consider spin-up electrons. 
Then, Eqs. (11) and (17) reduce to the ansatz 

MQ=ViMZ/^) ^ fft(0 = «3^n-iK/V2), 

(Al) 

valid for n > 1. For convenience, we have expressed the 
parabolic cylindrical functions D n (£) by the eigenfunc- 
tions of the one-dimensional harmonic oscillator, 

MO = D n (V2Z')/^riV^ = e- {a / 2 ff„(0//2»n!^, 

(A2) 

where H n (£') is the n-th Hermite polynomial. Inserting 
Eq. (Al) into Eq. (8) and using the recurrence relations 
for the parabolic cylindrical functions (9) and (10) leads 
to the eigenvalue problem 



C+M 



(T>+B)(2n+1) 
j? 

V2nA 
Ib 



2 



C-M 



V2nA 
Is 

(D-B)(2ra-1) 
p 



2 



= E 



(A3) 



By determining the eigenvalues of Eq. (A3) and their cor- and their respective (normalized) eigenstates 
responding eigenvectors, we find the Landau levels (25) 
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VI 



/ (V2^.A/iB=FAt,n/2)-[jM-(2Bn+P)/;^+(g.-g h )MBB/4] ^ 

y / A ti „(A ti „=F2v / 2n-4/if3)v / lB" 
(V2^A/l B T^, n /2) + [M-(2Bn+V)/l 2 B +{g a - gil )p, B B/4] M 



\/ A t,n(A t 



,q=2 v / 2H-4/iB)%/TB 






(y-ki 



(A4) 




2£n + £> 



9h 



(A5) 

Whereas Eqs. (Al) and (A4) are valid for n > 1, one can 
also choose n = to satisfy the boundary conditions. 
Instead of Eq. (Al), one then has the ansatz 



/t(O=«i0o(£/\/2) and . 9T (£)=0, 



(A6) 



which yields the single Landau level given by Eq. (26) 
and its corresponding (normalized) eigenstates 



*0,fc(^») 



/I 




Vo 



(A7) 



If a similar procedure is applied for the spin-down 
states, one finds the Landau levels given by Eq. (27) with 
the eigenstates 



Akx 



VI 



-(V2^A/l B ±A Un /2)-[M-(2Bn-V)/ll-(g c -g h ) l i B B/4] ( y-kij 

I — y^7i — X l 7 

v /A^,„(A ii „±2 v / 2Syl/iB)vT5" V 

-(v / 277^/iB±A 4 .,„/2) + [A^-(2Bn-P)/; B -( g o- fl h)A'B.B/4] , / y-fci 



y/ ' A ltn {^ l<n ±2^A/l B )^ 



(A8) 



where 



9 /2aA 2 / 2ffn-2? gc -g h 
Aj.,„ = 2\/ — V [M ^ — )J> B B 



f 2 



(A9) 

and the single Landau level given by Eq. (28) with the 
eigenstate 



e lkx 1 i /y-itn 



VI VU V i 




(A10) 



Appendix B: Ground-state magnetization 

As mentioned in Sec. IV A, the ground-state en- 
ergy (50) can be split in a — possibly not continuously 



diffcrcntiable — contribution from the uppermost valence 
band, fi dia (-B) given by Eq. (51), and a contribution from 
the remaining valence bands, Q,q(B) given by Eq. (52). 
Likewise, one can divide the magnetization of the ground 
state into 



M diB {B) 



i an di .(g) 

S dB 



and 



- 1 dQojB) 

Mo{B) --S^B— 



(Bl) 



(B2) 



Figure 23 shows the contribution to the magnetization 
from the uppermost valence band, Mdi s {B) , for parame- 
ters corresponding to the quantum-well thickness of d = 
7.0 nm, that is, the inverted regime. Here, one can clearly 
see the discontinuity of Mdi s (B) at B = B c . Compar- 
ing Mdis (B) to the non- vacuum contribution M (T, fx, B) 
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FIG. 23: (Color online) Magnetic field dependence of the mag- 
netization Mdi s (B) (corresponding to a quantum-well thick- 
ness of d = 7.0 nm). 



E =-150meV 




0.5 2 4 6 8 10 

B[T] 



FIG. 25: (Color online) Magnetic field dependence of the vac- 
uum magnetization Mo (B) corresponding to a quantum-well 
thickness of d = 7.0 nm and a = 10. 
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FIG. 24: (Color online) Magnetic field dependence of the 
non-vacuum magnetization M (T, n, B) (corresponding to a 
quantum-well thickness of d — 7.0 nm) for T = 10 K and 
different densities (n d = 10 14 , 10 15 , 10 16 1/m 2 ). 
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FIG. 26: (Color online) Magnetic field dependence of the vac- 
uum magnetization Mo (B) corresponding to a quantum-well 
thickness of d = 5.5 nm and a = 10. 



which is shown in Fig. 24 for T = 10 K and different den- 
sities illustrates how the discontinuity of M (T, fi, B) is 
canceled by the discontinuity of Mdi S (B) . The resulting 
magnetization can be seen in Fig. 20 in Sec. IV B. 

Apart from the contribution of M^s (B)+M (T, //, B), 
there is also a contribution arising from the remaining va- 
lence bands, Mq(B). When using the effective model for 
HgTe quantum wells given by Eq. (2), the valence band 
Landau levels are not bounded from below and, thus, the 
sum over them is divergent. However, the effective model 
used in this manuscript is only valid for low energies and 
there should be a lower bound for the valence band Lan- 
dau levels of the real band structure. To remedy this, 
we adopt the approach from Refs. 45-47 and introduce a 
smooth cutoff function g co (t) = E" D /(e a + E" Q ) which we 
include in the thermodynamical quantities to smoothly 
cut off the respective summation over the Landau levels. 
Here, E co and a denote the energy cutoff for the valence 
band Landau levels and a positive integer, respectively. 



Figures 25 and 26 show the contribution from Mo (B) for 
a = 10, several different energy cutoffs E co , and band 
parameters in the inverted (d = 7.0 nm) and normal 
(d = 5.5 nm) regimes, respectively. The main feature 
in these graphs is the decay of the magnetization with 
increasing magnetic field, indicating a negative suscepti- 
bility and therefore diamagnctism. 



Appendix C: Magnetization: Simplified model 

In the following, we briefly discuss the magnetization 
for the special case of the reduced model for Eq. (2) men- 
tioned in Sec. II. If one chooses C = 0, the bulk Landau 
levels (25)- (28) reduce to 

EVI(0) = ±M (CI) 
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and the degenerate levels 



(C2) 



in this case. 

If the simplified expressions (CI) and (C2) are 
used, the different contributions to the grand potential, 
a (T,n,B), Eqs. (48), (49), and (50), read as 



_. OO 

a (T,n,B) = -f(0) + jr / f(n), 



n=l 



and 

where 

/(*) 
and 



a (T, / i )J B) = n e (t,-h,b), 



I OO 

n (B) = -g(0) + J29(n), 



2SB 



In 



x e -p(y/2xA2/l 2 B +M*-v) 



2SB 2xA 2 ... 
g(x) = — — W-r*— +M 2 



/-' 

' r.-> 



(C3) 



(C4) 



(C5) 



(C6) 



(C7) 



In the following, we will look at the behavior of the 
magnetization in the regime of 2[Af3/l B ) 2 < 1 as well 
as the de Haas- van Alphen oscillations within the model 
given by Eqs. (CI) and (C2). For both cases, we assume 
to be in the degenerate limit, that is, 1. Since 

the Landau levels of this reduced model correspond to 
those of two-dimensional Dirac fermions, most notably 
those of (monolayer) graphene, one can apply the same 
procedures as in these cases. 



a. 'Weak' magnetic fields 

For magnetic fields with 2(A/3/l B ) 2 <C 1, we follow the 
classic Landau approach 48 and use the Euler-Maclaurin 
formula to express a (T, fi, B) as 



OO 

a {T,n,B)n J dx.f[x) 



1 df(x) 



12 dx 



(C8) 



When conducting the transformation x/l 2 B — > x, one can 
see that the integral in Eq. (C8) [denoted as F[T, n) in 
the following] does not depend on the magnetic field and 
one arrives at 



a (T,fi,B)^F(T,fx)- 



SA 2 



Y1-kI% \M\ c P(\m\-h) + i 



(C9) 



By the same procedure [and assuming a cutoff for g(x) 
we obtain 



a (B) « c 



SA 2 



\2nl% \M\' 



(CIO) 



where cq does not depend on the magnetic field. Then, 
the grand canonical potential can be written as 

a [T, (jl, B) = Q (B) + a (T, fjt, B) + a (T, -h, B) 
SnA 2 B 2 sinh(/3|X|) 



3$§ \M\ cosh(/3|7W|) +cosh(/3/^) ; 



(Gil) 

where the different _B-indepcndent contributions have 
been combined in the function F(T,n). Note, that 
the expansion used to arrive at Eq. (Cll) is valid for 
2{AP/l B ) 2 « 1. 

Consequently, we find for the magnetic susceptibility 



Xtot (T, n) 



2irA 2 



sinh(/9|M|) 



3$§ \M\ cosh [fi \M\) + cosh [fin) ' 

(C12) 

implying that the system is diamagnetic. This result gen- 
eralizes the zero-temperature formula of graphene found 
in Ref. 46, but also the M = model of Pbi_a;Sn x Te 
interface states found in Ref. 23. 



6. De Haas-van Alphen oscillations 

To calculate the de Haas-van Alphen oscillations for 
\fx\ > wc only need to look at the non- vacuum con- 
tributions a (T, fi, B) and a {T, fx, B). We again follow 
Ref. 48 as well as Ref. 49 and use Poisson's summation 
formula to write 



a(T,/z,B)« / dx/(a;) + 2Rc 



OO 

OO „ 

£ / dxf(x)c 2 * ikx 

fc=i i 



(C13) 

where the first and second terms describe the non- 
oscillating and oscillating parts of the grand potential, 
respectively. Here, we are interested in the oscillating 
part [denoted by 0° sc (T, /j,, B) in the following] . This 
part can be rewritten as 



a ac (T, (J,, B 
ASB 



k=l 



2nik£; 



\M/p\ 



where 



and ^ 



e [y-8gn(fi)]/£ _|_ I 

(C14) 



A 



y 



M 2 



(C15) 



i/GSH). 

We first consider the case \i > \M.\- In this 
major contribution to the integral originates from the 
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vicinity of the Fermi level, that is, from y ~ 1, whereas 
the integrand is damped for values y > 1. Therefore, we 
expand x(y) around y = 1 and replace the lower bound- 
ary of the integral by y — > — oo. Changing the integration 
variable to x = (y — l)/£, we find that the oscillating part 
of the grand potential is given by 



2SB 



He 



1C 



ikilB/Af^-M 2 ) 



x / d.T 



p 2iri{p,l B /A) £kx 



(C16) 



+ 1 



with fi > \A4\. For fi < \A4\, the contribution from 
the oscillating part of the electrons is much smaller than 
Eq. (C17) and in the case of fj, < — \M\, the main 
contribution arises from the hole contribution given by 
Q* ac (T, fj,, B) = Q° so (T, -fj,, B). Thus, the total oscillat- 
ing part of the grand potential is given by Eq. (CI 7) for 
any > By taking the derivative, one obtains 

the oscillating part of the total magnetization, which is 
periodic in 1/B. 



Computing the above integral, we can write the oscil- 
lating part of the electronic contribution to the grand 
potential as 



2SB ^ cos [irk{l B /A) 2 (m 2 - M 2 )} 



$ /3 



E 

k=l 



fcsinh [2n 2 k^(fil B /Ay 



(C17) 



Finally, we emphasize that this reduced model dis- 
cussed here cannot describe a transition between inverted 
and normal band structures and can thus only be used 
for magnetic fields well below the crossover point (or for 
situations where there is no crossover at all). 
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